Supersolid droplet crystal in a dipole-blockaded gas 
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A novel supersolid phase is predicted for an ensemble of Rydberg atoms in the dipole-blockade 
regime, interacting via a repulsive dipolar potential "softened" at short distances. Using exact 
numerical techniques, we study the low temperature phase diagram of this system, and observe an 
intriguing phase consisting of a crystal of mesoscopic superfluid droplets. At low temperature, phase 
coherence throughout the whole system, and the ensuing bulk superfluidity, are established through 
tunnelling of identical particles between neighbouring droplets. 

PACS numbers: 67.80.K-, 32.80.Rm, 67.85.Hj, 67.85.Jk, 67.85.-d, 02.70.Ss 



The search for novel phases of matter drives much of 
the current research in condensed matter physics. Of 
particular interest are phases simultaneously displaying 
different types of order. A chief example, of great cur- 
rent interest, is the so-called supersolid^ namely a phase 
featuring crystalline order, and also capable of sustaining 
dissipation-less flow. Attempts to observe experimentally 
a supersolid phase of matter, primarily in a crystal of 
solid helium, have spanned four decades since early the- 
oretical predictions [1] . The most credible claim of such 
an observation to date [H |3], has been subjected to in- 
depth scrutiny over the past few years, and it seems fair 
to state that agreement is lacking at the present time, as 
to whether experimental findings indeed signal a super- 
solid phenomenon [4 . 

A new, fascinating avenue to the observation of super- 
solid and other phases of matter not yet observed (or 
even thought of), is now opened by advances in cold 
atom physics, providing not only remarkably clean and 
controlled experimental systems, but also allowing one 
to "fashion" artificial inter-particle potentials, not aris- 
ing in any known condensed matter system. This allows 
one to address a key theoretical question, namely which 
two-body interaction potential(s), if any, can lead to the 
occurrence of a supersolid phase in free space (i.e., not 
on a lattice). 

In a recent article [5], Henkel et al. have proposed, 
based on a mean- field treatment, that a Bose conden- 
sate of particles interacting through an effective poten- 
tial which flattens off at short distance, might support 
a density modulation. In this Letter, we show by first 
principle numerical simulations that interaction poten- 
tials which combine a long-distance repulsion with a 
short-distance cutoff, lead in fact to the appearance of a 
novel self-assembled crystalline phase of mesoscopic su- 
perfluid droplets in a system of bosons. Furthermore, 
such a crystal can turn supersolid in the T ^ limit, 
as tunneling of particles across neighbouring droplets 
takes place, and superfluid phase coherence is established 
across the whole system, as individual separate Bose con- 



densates (droplets) organize into a single, global conden- 
sate. Specifically, we consider the following two-body 
potential: 



v{r) 



D/a^ if r < a 
D/r^ if r > a 



(1) 



D being the characteristic strength of the interaction. 
This kind of interaction potential can be realized with 
cold dipole-blockaded Rydberg atoms [6UTQ]. The pa- 
rameters D and a above can be controlled with external 
fields [5l[TTj (we come back to this point below). 

Our system of interest comprises N identical bosons of 
mass m^ confined to two dimensions [21 . The many-body 
Hamiltonian is the following (in dimensionless form): 



n 



1 ^ 






(2) 



where r^j = Ir^ — rj| is the distance between particles i 
and j, and v is given by Eq. ([l]). All lengths are expressed 
in terms of the characteristic length Tq = mD/h^^ and we 
introduce a dimensionless cutoff Re = a/vo for the poten- 
tial (IT]). The system is enclosed in a square cell of area A, 
with periodic boundary conditions. The particle density 
is n = N/A, but we shall express our results in terms of 
the (dimensionless) inter-particle distance rs=l/\/nr^. 
The energy scale is eo = D/r^ = h^ /mrl. 

The low-temperature phase diagram of such a system 
has been explored by means of first principles numerical 
simulations, based on the Continuous-space Worm Al- 
gorithm [12]|13]. Numerical results shown here pertain 
to simulations with a number of particles N varying be- 
tween 50 and 400, in order to carry out extrapolation 
of the results to the thermodynamic limit. Our ground 
state estimates are obtained as extrapolations of results 
at finite temperature. Details of the simulations are stan- 
dard, as the use of the potential ([l]) entails no particular 
technical difficulty. 

In the limit Re <C r^, the truncation of the dipolar 




FIG. 1: Snapshots of a system of bosons interacting via po- 
tential (IT]), at the four different temperatures 200 (a), 20 (b), 
1.0 (c) and 0.1 (d), expressed in units of Co- Points shown are 
taken along individual particle world lines. The nominal value 
of rs in this case is 0.14, whereas the cutoff of the potential 
Q is Rc=0.3. 
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potential at short distances does not play an important 
role, and the low temperature phase diagram of ^ is 
that of purely dipolar bosons in two dimensions, investi- 
gated previously by several authors [14[ [15] . It is known 
that for Tg < r^ = 0.06 the ground state of the system 
is a triangular crystal, whereas for r^ > rf" = 0.08 it is 
a uniform superfluid (in the intermediate density range 
a more complex scenario is predicted [H]). As we show 
below, a very different physics sets in when Re > r^, in 
the density ranges which correspond to either the crys- 
talline or superfluid phase in the purely dipolar system. 
Fig. fl] shows typical configurations (i.e., particle world 
lines) produced by Monte Carlo simulations of a system 
of bosons interacting via the potential ([l]) , at a nominal 
density corresponding to r^ = 0.14, at different temper- 
atures spanning three orders of magnitude. The value of 
the cutoff Re in this case is 0.3. At the highest temper- 
ature, a simple classical gas phase is observed, as shown 
by the pair correlation function ^(r), shown in Fig. ^(a), 
which is just a constant (note that g{r) does not vanish 
at the origin, owing to the flattening off of the potential 
at short distance). As T is decreased, an intriguing ef- 
fect takes place, namely particles bunch into mesoscopic 
droplets, in turn forming a regular (triangular) crystal. 
This is shown qualitatively in the snapshots in Fig. [l] but 
also confirmed quantitatively by the structure of the g{r) 
as well (Fig. |2]^a)), which displays pronounced, broad 
maxima, as well as well-defined minima, where the func- 
tion approaches zero. We henceforth refer to this phase 
as the droplet- crystal phase. 



FIG. 2: (Color online) Results shown are for rs=0.14 and 
i?c=0.3. Temperature is in units of Co- (a): Pair correla- 
tion function g(r) at a temperature T= 200 (triangles), 20 
(squares), 1.0 (diamonds) and 0.1 (circles). The simulated 
system comprises N=200 particles, (b): Superfluid density 
vs. T for systems with A/'=100 (square), and 200 (diamond) 
particles, (c): Frequency of occurrence of permutation cycles 
of length L at the same four temperatures reported in panel 
(a). Longer permutation cycles occur at lower temperature. 



The formation of such droplets is a purely classical 
effect, that depends on the flattening off of the repul- 
sive inter-particle potential below the cutoff distance. In 
fact, a simple estimate of the number N^ of particles per 
droplet, can be obtained by considering a triangular lat- 
tice of point-like dipoles, each one of strength ex N^ (as 
it comprises N^ particles), and by minimizing with re- 
spect to Nd the potential energy per particle, for a fixed 
density. The result is 



N. 






(3) 



where 7 ~ 2.79. Eq. ^ furnishes a fairly accurate esti- 
mate of Nd for the (wide) range of values of the param- 
eters rs and Re explored here. For instance, using the 
parameters of Fig. [l] we find from ([3| N^ ^ 13, which 
agrees quite well with our simulation result. It is worth 
noting that a similar sort of pattern formation, due to 
competing interactions, has been previously established 
for classical colloidal systems [171 UHl • 



IDC C SF 

I ■ Re < rs^ 



FIG. 3: Schematic ground state phase diagram of (l2| as a 
function of Ts- The superfluid droplet crystal (SDC) is sand- 
wiched between an insulating droplet crystal (IDC) and a su- 
perfluid (SF). For Re '^ T5 , an IDC, a single-particle crystal 
(C) and a superfluid phases are observed. The widths of the 
SDC and C regions depend on the value of Re. 



In the T ^ limit, long exchanges of identical parti- 
cles can take place, as a result of particles tunneling from 
one droplet to an adjacent one. Long exchanges of parti- 
cles can result in a finite superfluid response throughout 
the whole system, and indeed for R^ > r^ we observe 
such a bulk superfluid signal, in a range of values of r^ 
in the vicinity of Rc/"^- A typical result is shown in Fig. 
[2|b) p2 . Because superfluidity arises in concomitance 
with the droplet-crystal structure, the denomination su- 
per solid seems appropriate. 

In order to establish that droplets are individually su- 
perfluid, one may consider the statistics of permutation 
cycles. Fig. El (c) shows the frequency of occurrence of 
exchange cycles involving a varying number L of parti- 
cles (1 < L < A^), at three different temperatures, at 
the physical conditions of Fig. [l] As one can see, as the 
temperature is lowered exchange cycles involving grow- 
ing numbers of particles occur. At low temperature they 
involve almost all the particles in the system; however, 
even at a higher temperature (e.g., T=20 in Fig. [2|c)) 
one observes exchanges comprising a number of particles 
up to ~ Nd^ i.e., particles inside an individual droplet. 
This is evidence that droplets are individually Bose con- 
densed and superfluid, even though the system as a whole 
does not display superfluidity. That droplets should be 
superfluid at low T is not surprising, given that particles 
in a droplet are essentially non-interacting, due to the 
flatness of the potential at short distance. However, that 
droplets are themselves superfluid does not imply that 
a bulk supersolid phase will always occur in the T ^ 
limit, as discussed above. 

At T = 0, the supersolid phase is sandwiched be- 
tween an insulating droplet crystal at high density (i.e., 
lower Vs) and a homogeneous superfluid phase at lower 
density. For Re '^ r^, only two insulating phases are 
observed, namely the insulating droplet crystal at high 
density and the crystal of single particles, already de- 
tected in Refs. [HI [15], as well as a superfluid phase at 
lower density. All of this is summarized in the schematic 
phase diagram shown in Fig. [3] It is important to stress 
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FIG. 4: (Color online) Left panels: Monte Carlo snapshots of 
a system oi N — 400 particles, interacting via the potential 
([1]), confined in a harmonic trap of strength F = 500 eo, at 
the two temperatures T=100 eo (a) and T=0.5 eo (b). Right 
panels: corresponding momentum distributions, all normal- 
ized to unity for comparison purposes. The value of Re in 
this case is 0.3. The development of secondary peaks at low 
temperature signals the occurrence of a supersolid phase. 



that supersolid behaviour in this system originates from 
tunnelling of particles between droplets which are them- 
selves individually superfluid, so that the individual su- 
perfluid droplets connect to form a bulk superfluid. This 
is reminiscent of the phase-locking mechanism in a (self- 
assembled) array of Josephson junctions. 

The results discussed so far pertain to numerical sim- 
ulation of the system described by Eq. Q in its bulk 
phase. However, in any experiment aimed at probing 
the physics of such a system, the assembly of particles 
must necessarily be finite (a few thousand particles is a 
typical number for current experiments with cold dipo- 
lar atoms), confined by an external potential. In order 
to enable a direct comparison with possible future ex- 
periments, we have performed simulations of the same 
system spatially confined in-plane by a harmonic trap, 
i.e., the term T^.rf is added to Eq. (|2|, F = muj^/2 
being the strength of the trap. 

Fig. [4] shows typical many-particle configurations of 
a trapped system comprising A^=400 particles, at two 
different temperatures. Also shown are the associated 
momentum distributions n(/c), which are obtained by 
Fourier transforming of the spherically and translation- 
ally averaged one-body density matrix, computed by 
Monte Carlo. 

Here too, droplets with a well-defined average number 
of particles form, and organize themselves on a triangu- 
lar lattice. Correspondingly, the momentum distribution. 



which is directly observable experimentally by time-of- 
flight measurements [T9], develops a sharp central peak, 
with additional structure on its sides. The secondary 
peaks correspond to oscillations in the one-body density 
matrix, in turn reflecting particle tunnelling to adjacent 
droplets. They are therefore connected to the appearance 
of the supersolid phase, as explained above. 

Summarizing, accurate numerical simulations of a sys- 
tem of dipolar particles interacting via a potential soft- 
ened at short distance, reveal the existence of a low tem- 
perature crystalline phase of superfluid droplets. This 
phase turns superfluid (supersolid) at T ^ through 
a mechanism of tunnelling of particles between adjacent 
droplets. At higher density this tunnelling is suppressed 
and an insulating droplet crystal occurs, a phase which 
has not previously been predicted. The interaction that 
underlies such intriguing, until now unobserved physi- 
cal behaviour, can be realized with dipolar atoms in the 
dipole-blockade regime. 

A comment is on order, concerning the dependence of 
the results on the particular form of potential utilized 
here, namely Eq. M with its abrupt, sharp cutoff at 
r = a. First off, the superfluid droplet crystal phase 
does not crucially depend on the dipolar form of the in- 
teraction at long distances. Indeed, it is also observed 
in our simulations for a van der Waals-like potentials 



(i.e., v{r) 



^). Secondly, we have obtained qualita- 



tively similar results with different model potentials, fea- 
turing a more realistic "flat" region at short distances, 
as well as a smoother merge of long- and short-range 
behaviours. For example, we considered the potential 
V{r) = I)/(a^ + r^), which is naturally realized in a 
cold gas of alkali atoms by weakly dressing the ground- 
state 1^) of each atom with an excited Rydberg state \r) 
with a large dipole moment d, in the kDebye range [11 . 
While several dressing schemes are possible [SI [20], here 
we consider \r) as the lowest-energy state of a Rydberg 
manifold with principal quantum number n for an atom 
in the presence of a homogeneous electric fleld F < Fit 
in the linear Stark regime, with Fit cx 1/n^ the Inglis- 
Teller limit. For a laser with (effective) Rabi frequency 
Q and red detuning |A| ^ Q^ the dressed groundstate 
reads \g) - \g) - {n/2A)\r), and thus D ex (l^/A)^^^, 
with d (X n^^ while the cutoff a ^ {cP /h\A\eo)^^^ arises 
because of the Rydberg-blockade mechanism [6 . Spon- 
taneous emission rates 7^ from \r) are strongly reduced 
to values of at most 7 c^ {Q/A)'^jr- Observability of 
the phases above beneflts from large values of eo, and 
in particular of the ratio eo/{hj) = [Q'^ /{47r\A\jr)]Rc^ 
favoring comparatively small values of To [23 . For exam- 
ple, for ^^Rb atoms in an electric fleld F = 25kV/m, 



with n = 20, ^/27r = 50MHz and |A|/27r = 3GHz, 
we obtain d = 1450 Debye, a c^ 400 nm, To ^ 210 
nm, Co/ikj) 0:^ 90, and Co/kB — 120 nK. Collective 
many-body effects in the Rydberg-blockade regime ^20] 
not described by Eq. ^ should be negligible provided 
{Rc/tsY <C (2A/r^)^, which is readily satisfled for pa- 
rameters as in Fig. [T] 
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Confinement to two dimensions can be achieved using a 

tight (magnetic or optical) trap along z. 

We compute directly the hulk superfluid fraction through 

the usual winding number estimator for the bulk system. 

2D scattering requires Vo > a±, with a± the transverse 

harmonic oscillator length, in the tens of nm range. 



